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CHANGE OF BASIS FOR GENERALIZED CHEBYSHEV 
KOORNWINDER’S TYPE POLYNOMIALS OF FIRST KIND 

MOHAMMAD A. ALQUDAH AND MAALEE N. ALMHEIDAT 


Abstract. In this paper an explicit form of generalized Chebyshev Koorn- 
winder’s type polynomial of first kind in terms of the Bernstein basis of fixed 
degree n is provided. Moreover, we investigate generalized Chebyshev Koorn- 
winder’s type polynomials of first kind and Bernstein polynomials change of 
bases. 


1. Introduction 


Approximation is important to many numerical methods, as many arbitrary 
continuous functions can be approximated by polynomials. On the other hand, 
polynomials can be characterized in many different bases. Every type of polynomial 
basis has its strength, where by suitable choice of the basis, numerous problems can 
be solved and many complications can be removed. 

Characterization of the generalized Chebyshev Koornwinder’s type polynomials 
of first and second kind using Bernestin basis were discussed in BIZ] respectively. 
On the other side, Rababah [7] investigated the transformation of classical Cheby¬ 
shev Bernstein polynomial basis. 

In this paper, we generalize the results in [7], where an explicit closed form 
of generalized Chebyshev-Koornwinder’s type polynomials of first kind (ChK-1) 
of degree r < n expressed in terms of Bernstein basis of hxed degree n, and a 
closed forms of generalized Chebyshev Koornwinder’s polynomials and Bernstein 
polynomials bases transformations are provided. 

1.1. Bernstein polynomials. The Bernstein polynomials have been studied thor¬ 
oughly, there exist many great enduring works on theses polynomials [4]. The 
n -I- 1 Bernstein polynomials BJ^{x) of degree n, x G [0,1], are defined by B'^{x) = 



Bernstein polynomials are the standard basis for the Bezier curves and surfaces 
in Computer Aided Geometric Design (CAGD). Many applications that involve two 
or more Bezier curves require all involved Bezier curves to have the same degree, so 
we use the degree elevation defined by [S] to write Bernstein polynomials of degree 
r where r < n in terms of the Bernstein polynomials of degree n 



( 1 . 1 ) 


Although higher degree Bezier curves require more time to process, they do have 
more flexibility for designing shapes, where the basis are known to be optimally 
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stable. So, it would be very helpful to increase the degree of a Bezier curve without 
changing its shape mw,- 

Farouki [1] surveyed Bernstein basis properties, he described some of the key 
properties and algorithms associated with the Bernstein basis. Bernstein polyno¬ 
mials remarkable properties [4] make them significant for the development of Bezier 
curves and surfaces in CAGD. However, the Bernstein polynomials are not orthog¬ 
onal and could not be used effectively in the least-squares approximation [8]. the 
calculations performed in finding the least-square approximation polynomial of de¬ 
gree m do not decrease the calculations to obtain the least-squares approximation 
polynomial of degree to -|- 1. Since then, the method of least square approximation 
accompanied by orthogonality polynomials has been presented and examined. 

In the following, we define least-square approximations of a continuous function 
f{x) by using polynomials with standard power basis l,x,x^,..., x". 

Definition 1.1. For a continuous function f{x) defined on [0,1], the least square 
approximation requires finding a least-squares polynomial p^i^) = 
that minimizes the error E{ao,ai,... ,an) = — p'^{x)]'^dx. 

A necessary condition for E{ao, oi,..., a„) to have a minimum over all values 

OE 

oo, ai,..., a„, is 0 = ——, k = 0, ... ,n. Thus, for z = 0, 1,..., n, Ui that minimize 
oak 

Wfi^) - ELo«fc<^fe(a^)ll2 satisfy 

pi n „i 

/ f{x)4>i{x)dx = '^ak 4)k{.x)(l)fix)dx, 
do do 

which gives a system of (n 1) equations in (n 1) unknowns: a^, z = 0, ..., n, 
called normal equations. Those (n -I- 1) unknowns of the least-squares polynomial 
Pnix), can be found by solving the normal equations. By choosing (j)i{x) = x\ a.s a 
basis, then 

„i n „i n 

/ f{x)x^dx = ^Ofe / F+'^dx = ■ , ?, 1 ■ 

•^0 fe=0 k =0 

The coefficients matrix of the normal equations is Hilbert matrix, Etij = (z -|- 
j — 1)“^, i,j = l,...,n, which has round-off error difficulties and notoriously 
ill-conditioned for even modest values of n. However, such computations can be 
made computationally effective by using orthogonal polynomials. Thus, choosing 
{4iq{x), 4>i(x ),..., (pnix)} to be orthogonal simplifies the least-squares approxima¬ 
tion problem. The coefficients matrix of the normal equations will be diagonal, 
which simplifies calculations and gives a closed form for ai,i = 0,l,...,n. More¬ 
over, once Pn{x) is known, it is only needed to compute a„+i to get pn+i{x), which 
turns out to be computationally efficient, to get pn+i{x), which turns out to be 
computationally efficient. See [5] for more details on the least squares approxima¬ 
tions. 


1.2. Factorials and double factorials. The double factorial of an integer to 
is given by {2m — 1)!! = (2to — 1)(2to — 3)(2to — 5)... (3)(1) if to is odd, and 
to!! = (to)(to — 2)(to — 4)... (4)(2) if to is even, where 0!! = (—1)!! = 1. 
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Definition 1.2. For an integer n, the double factorial is defined as 


( 1 . 2 ) 



if n is even 
if n is odd 


From the definition o and the fact (1/2)! = we can derive the factorial 

of an integer minus half as 



n\(2n — 1)1! 
(2/^)!! 


1.3. The generalized Chebyshev-I Koornwinder’s type polynomials. For 

M, IV > 0, the generalized Chebyshev-I Koornwinder’s type polynomials [1] of 
degree n 


(1.3) 




(2n- 1)1! 
(2n)!! 


Tnix) + ^ 


miXk 

22 '=( fc !)2 


Tk{x), 


n = 0,..., 00 , 


^ 

are orthogonal on the interval [—1,1] with respect to the measure - dx + 

TT 

Md-i + N6i, where 6x is a singular Dirac measure, 

Tn{x) are the Chebyshev polynomials of the first kind (Chebyshev-I), and 


(1.4) 


Xk — 


4M 


41V 


AMN 


{2k-i){k-l)\ (2/fc-3)(fc-1)1 (fc-l)!(fc-2)!’ 


The following theorem [T] shows how generalized Chebyshev-I Koornwinder’s 
type polynomials X7r’^'^\x) of degree r can be written explicitly as a linear com¬ 
bination of the Bernstein polynomial basis of degree r. 


Theorem 1.1. [I]/M. AlQudah) For M,N > 0, generalized Chebyshev-I Koorn¬ 
winder’s type polynomials XFr^'^\x) of degree r have the following Bernstein rep¬ 
resentation, 


(1.5) 


(2r- 1)!! 
(2r)!! 


^(-l)^-*i7..,i3[(x) + ^ 

i=0 k—0 


m\Xk 

22 fc ( fc !)2 


k 




where Xk defined in (HD, V^,r 
Moreover, the coefficients rji^r 


(2r)(2,-) 

= C = 0A,---,r, and Po^r 

satisfy the recurrence relation 


( 1 . 6 ) 


'^i,r 


{2r-2i + l) 

(2*-l) 


Vi-ljr! 


i = l,.. 


r. 



2. Change of Bases 

In this section we provide a closed form for the transformation matrix of gen¬ 
eralized Chebyshev-I Koornwinder’s polynomials basis into Bernstein polynomials 
basis, and for Bernstein polynomials basis into generalized Chebyshev-I Koorn¬ 
winder’s polynomials basis. 
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2.1. Generalized Chebyshev-I to Bernstein. Rababah [7] provided some re¬ 
sults concerning the classical Chebyshev case. In the following theorem we general¬ 
ize the procedure in [7] . The following theorem will be used to combine the superior 
performance of the least-squares of the generalized Chebyshev-I polynomials with 
the geometric insights of the Bernstein polynomials basis. 


Theorem 2.1. The entries = 0,1,...,n of the matrix transformation 

of the generalized Chebyshev-I polynomial basis into Bernstein polynomial basis of 
degree n are given by 
( 2 . 1 ) 


Mf 


(ry^ (2r-)! 
\i) 22’'(r!)2 


min(i,r) 

E 

/c=max(0,2+'r—n) 







-1 


(2fc)! 

22fe(fc!)- 


min(2,/c) 

■ E 


(- 1 ) 


k-j 


n — k 

i- j 


k-\ 



Proof. Write Pn{x) as a linear combination of the Bernstein polynomial basis as 
and as a linear combination of the generalized Chebyshev-I polynomials as follows 
Pn{x) = TTr^.CrB-ix) = d.^^'^\x). 

We are interested in matrix M, which transforms the generalized Chebyshev-I 
coefficients Bernstein coefficients Ci = X)r=o 

By using Theorem 1 1.1 1 to write the generalized Chebyshev-I polynomials (11.31) of 
degree r < n in terms of Bernstein polynomial basis of degree n as 


( 2 . 2 ) 


= E ^ = 0,1, 


i=0 


where the entries of the (n -I- 1) x (n -|- 1) basis conversion matrix N. Thus, the 
elements of c can be written in the form 

C" — \ rj 

Cl — Z^r=0 

Now, we need to write each Bernstein polynomial of degree r where r < n in 
terms of Bernstein polynomials of degree n. By substituting the degree elevation 
dni) into dni) and rearrange the order of summations, we find that the entries of 
the matrix N = M'^ for r = 0,1,..., n are given by 


AT". = 


(2r)! 


min(i,r) 

22 r ( j , n 2 E (“1) 

fe=max(0,2+r—n) 


r — k 


n — r\ (r — 
i — k 


' 2 
r — k 






(2fc)! 


m.in(i,fc) 

22 fe ( fc !)2 E (“1)^ ^ 

j=max(0,2+A;—n) 


n-k^ y-\ 
3 



□ 

In the following Corollary, we express the generalized Chebyshev-I polynomials 
of degree r < n in terms of the Bernstein basis of fixed degree n. 
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Corollary 2.1. The generalized Chebyshev-Ipolynomials ..., 

of degree less than or equal to n can be expressed in the Bernstein basis of fixed de¬ 
gree n by the following formula 

n 

=^NTBf'{x), r = 0, 

i=0 

where = /x"^ + X]fc=o ^klJ-k,i Tr,i defined as 


(2r)! 


min(i,r) 

b''r,i = 22r('^n2 X/ 

fe=max(0,i+r—n) 


(-ly 


/n—r\ /2r\ /2r\ 
fc \i—k) \ r ) \2k) 

2“(”) ' 


Proof. From equation (12.21) From the proof of Theorem l2.ll any generalized Chebyshev- 
I polynomials J^r^’^\x) of degree r < n can be expressed in terms of Bernstein 
basis of fixed degree n as t^r^’^\x) = ^^^=0 N"^B^{x), r — 0 , 1 ,..., n, where 
the matrix N can be obtained by transposing the entries of the matrix M defined 
in (12.11) . But, 


' 2 
r — k 


(2r- 1)!! 


(2r- 1)!! 


-yr — k 


2'-(r - fc)! (2/c - 1)!! 2^k\ (2r-2fc-l)!! 

1 (2r-l)!! (2r-l)!! 


2'-(r - k)lkl (2k - 1)!! (2(r - k) - 1)!!' 
Using the fact (2r)! = (2r — l)!!2’'r! we get 

/„_rWr-iwr-U /n-rW2rW2r\ 

\i — k/ \r—k/ \ k / \i—k) \ r / \2k/ 


(2.3) 


C) 


22.(«) 


So, the entries can be rewritten as 

r 


k^O 


where 


(2r)! 


min(i,r) 

l^r.i = 22r('^n2 X/ 

fe=max(0,z+r—n) 


(- 1 ) 


r — k 


/n—r\ (2r\ /2r\ 
\i-k) \r ) \2k) 

2“(”) 


□ 


2.2. Bernstein to generalized Chebyshev-I. Farouki [3] discussed some ana¬ 
lytic and geometric properties for Bernstein polynomials. It is worth mentioning 
the product of two Bernstein polynomials is a Bernstein polynomial and given by 

Bff(x)Bffi(x) = (x), and the Bernstein polynomials can be integrated 

(, t+j 1 ■' 


easily as fg BJt(x)dx = k = 0,1,... ,n. In addition, we have the following 

Lemma that will be used in the proof of next theorem. 


Lemma 2.1. [I](^M. AlQudah) Let he the Bernstein polynomial of degree 

n and (x) be the generalized Chebyshev-I polynomial of degree i, then for 
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i,r = 0,1,... ,n we have 


P ^ 

/ {1 — 
Jo 


n\ {2i)l 


22i(i!)' 


E(-i) 


i—k 




*- 2 
i — k 


1 


1 , 


B{r + k+ -,n + i — r — k+-) 


+ 


±N(:)5^g(-ir 


d=0 


d — ^d — 


d-jj 


]B{r+j + -,n + d-r-j + -) 


where B{x,y) is the Beta function. 


Proof. See [1] for the proof. 


□ 


Theorem 2.2. The entries , i, r = 0,1,..., n of the matrix of transformation 

of the Bernstein polynomial basis into the generalized Chebyshev-I polynomial basis 
of degree n are given by 

2^ , /^ ^ (2rf)!A, ; 

(2z)! V (2^)! J 22'i(d!)2 ’ 


where Xi defined in (ITrl) . = 


r 2/tt if i = 0 
\ I/tt ifij^O ’ 


(2.4) = 




(1 + A02 Vr 


*;:r = i : 


i—k 




22* 


B{r + fc+i,n + i- r- A: + i), 


and B(x,y) is the Beta function. 


Proof. To write the Bernstein polynomial basis into generalized Chebyshev-I poly¬ 
nomial basis of degree n, we invert the transformation Ci = X]"=o Mf.^dr to get 

(2.5) d = M-^c. 

Let ^\ r = 0,..., n be the entries of M~^ and N~^ respectively. The 
transformation of Bernstein polynomial into generalized Chebyshev-I polynomial 
basis of degree n can then be written as 

n 

( 2 . 6 ) Bf{x)=J2NT^^^^’^\x). 

i=0 


The entries , i,r = 0,1,...,n can be found by we multiply (12.611 by (a: — 

x'^)~^ and integrate over [0,1] to get 

(2.7) 

n 1 ^ ^ 1 

/ {x-x^)-"^B'^{x)3^^^'^\x)dx = ^Nft~P / {x-x^)-"^Sr}^'^\x)3r^^'^\x)d: 

Jo ’ Jo 


2=0 


By using the orthogonality relation 



x) ^x ^T„(x)Tm(x)dx 


0 if TO n 

^ if TO = n = 0 

TT if TO = n = 1, 2,... 


we get 
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/ (x — x^) ^ B^(x)^/^’^\x)dx = 

Jo 

Taking into account Theorem 12.11 we obtain 

$*(”) ^22*(i!)2 ^ (-l)*-fe f2i\ (2i 

2k 



( (2i)! ^ 


\ ( 

^ ( 2 i)! y 



TV" / = 


(l + AO^^ (2z)! 22* 


■E 


B{r -\- k + —,n + i — r — k + —) 


22*(*!)* 

( 2 *)! 


^ (2d)!Ad A (-1)'^”^ (‘Jd\ f2d 

0-2d(M\2 2^ 


22dM!)2 ^ 22*^ 

d=0 ^ ' j=0 


" )S(r + j + i,n + d-r-j + i)]. 


By reordering the terms we have 


(1 + A,)2(r) 

2/7r if i = 0 


22*(d)^ 


where = 


I/tt if j ^ 0 


-+ 

(2i)! 

, Ai defined in (11.411 


2j 

22*(d)2^' * 


( 2 *)! 


(2d)!Arf 


d=0 


22rf(d!) 


2 j,d 


^ (-1)^- ^ f2^\f2^ 

^k,i / ^ 2‘2'i 


2k 


B{r -\- k + -,n + i — r — k + -), 


and B{x, y) is the Beta function. The entries of M ^ are obtained by transposition 


of TV 


-1 


□ 
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